Connecting diffusion and dynamical heterogeneities in actively deformed amorphous 
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In this letter we explore the relations between tracer diffusion and flow heterogeneities in flowing 
amorphous materials. On the basis of scaling arguments as well as an extensive numerical study 
of an athermal elasto-plastic model, we show that there is a direct link between the self-diffusion 
coefficient and the size of cooperative regions at low strain rates. Both depend strongly on rate 
and on system size. A measure of the mean square displacement of passive tracers in deformed 
amorphous media thus gives information about the microscopic rheology, such as the geometry of 
the cooperative regions and their scaling with strain rate and system size. 
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Under the application of sufficiently large stresses, 
solid amorphous materials, such as dense emulsions 
colloidal and granular systems Q or molecular glasses 
3, undergo plastic deformation and flow. This yielding 
behavior usually comes along with peculiar spatial fea- 
tures, such as deformation heterogeneities, shear bands 
[3| or nonlocal flow behavior [l[. At the fundamental 
level, the emerging picture to describe these properties 
involves local dissipative rearrangements, the so called 
'shear transformations' introduced by Argon cascad- 
ing via long range elastic interactions. This scenario has 
been largely confirmed by a number of numerical 0, 0] 
and experimental Q studies of quasistatic deformation in 
amorphous systems. It was also fruitful on the theoret- 
ical side, and served as the basis for various rheological 
descriptions, which aim at describing the statistical col- 
lective properties of plastic events p, [l3|- Yielding, 
nonlocality in the flow behavior, enhanced mobility |3| 
and heterogeneities are therefore only different aspects 
of the same underlying physics. 

In this letter, we consider the question of tracer diffu- 
sion in such flowing materials. On the experimental side, 
results indicate that the diffusion of tracer particles, e.g. 
in sheared foam systems and in colloidal glasses, is strain 
rate dependent ll,l3. Further observations concern the 
strong nonlocal effects in the mechanical response even in 
quiescent regions far away from the flowing material [l^ . 
The observation of a coupling between enhanced particle 
diffusion and an imposed deformation is quite natural in 
the context of slowly driven systems, where the exter- 
nal driving provides a source of energy that may share 



some characteristics of a thermal bath [IJ]. Here, we 
investigate this connection based on the specific defor- 
mation mechanisms associated with the flow behavior. 
Tracer diffusion is expected to take its origin in the me- 
chanical noise induced by the spatio-temporal collective 
behavior of local plastic events. In this context, we first 



establish a direct link between flow induced dynamical 
heterogeneities and the diffusive dynamics of tracer par- 
ticles immersed in a deformed system (see Fig. [T|). We 
summarize this in the form of a general scaling relation 
between the diffusion coefficient and the 4-point corre- 
lations that quantify flow heterogeneities. This allows 
us to rationalize the strain rate dependence, as well as a 
strong flnite-size dependence of the diffusion coefficient in 
sheared amorphous materials. We first discuss the gen- 
eral scaling arguments for such a connection, and in a 
second part we support these arguments on the basis of 
large scale simulations of an elasto-plastic model of yield 
stress materials. 




FIG. 1: (color online) Schematic view on the heterogeneously 
distributed cascading plastic events (on the left) and the re- 
sulting long range effects on the diffusion of a tracer particle 
(on the right). 

A scaling argument. - We start from the microscopic 
picture of 'shear transformation events'. Each event is 
assumed to correspond to a local Eshelby transformation 
15|, with a strain amplitude Aeq localized in a zone of 
linear size a. The displacement field for an event oc- 
curring at the origin decays algebraically; in 2D it reads 
u(r) = (2a^Aeoxy/7rr'^) x r. If we consider the non-affine 
part of the dynamics of tracer particles in an athermal 
medium, their relative position tt will be determined by 
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the time integrated velocity field v caused by the sur- 
rounding plastic events, 



with r^-* = r{t = 0), the position of the ith plastic 
event and Np(t) the number of plastic events at time t. 

To give an estimate for the mean square displacement 
of these tracer particles in the highly cooperative regime, 
we introduce the notion of cooperative events that span 
a cooperative volume 14, and have a typical duration 
t* . The typical mean square displacement due to an 
cooperative event (avalanche, slip-line or more complex 
object, depending on the geometry) can be estimated 
as (A/2) p2 ^dj,j^ (fir'C{Y - r'), with p the den- 
sity of elementary events in the cooperative volume and 
C(r - r') = (u(r)u(r')). This yields 



N^(t') 
dt' ^ 



1=0 



v(rT(0-r,(i')) (1) 



(2) 



up to logarithmic corrections originating from the dis- 
placement field correlations The number of the 
cooperative events, N^, occurring during a strain inter- 
val A7 can be estimated by the total number of plastic 
events Np{A'-f) = L'^ A-f / Aeq divided by the number 
of events in the cooperative volume, Nc{A'y) = Np/pVc- 
If we further assume that the cooperative events are sta- 
tistically independent, the diffusion coefficient D for the 
tracer particles in the limit of weak strain rate 7 can be 
written as D = 7!) with 



D = —n. ~ — pVc . 
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The appropriate tools for describing cooperativity are the 
'4-point' correlation functions, that quantify the corre- 
lation in space and time between events, and are de- 
scribed in more detail below. The 4-point function of 
interest here will be the one associated with plastic ac- 
tivity, rather than particle displacements. The integral 
over space of such correlation functions, denoted by Xi{t) 
goes in general trough a maximum at a finite time i*, 
which can be used as a measure of the volume over which 
events are correlated. Therefore we expect, in macro- 
scopic systems that the diffusion coefficient should scale 
in the same manner, i.e. 



(4) 



As we will demonstrate in the example below X'^{^*) is not 
only strain rate but also strongly system size dependent. 
This behavior originates from the long range character of 
the displacement field created by the plastic events. As a 
result we recover the same strong nonlocal effects for the 
diffusion coefficient. This is consistent with the experi- 
mental observations, that predict a sensitive mechanical 



response to distant active regions [l3|, and with results 
from simulations of simple particle systems [l6| . 

An elasto-plastic mesoscale model - To test the ideas 
described above, we performed large scale simulations 
17| of a simple mesoscopic model for the flow of yield 
stress materials, introduced by Picard et. al. [l^. The 
model is restricted to athermal systems and to a plane 
geometry. It features the dissipation due to local plas- 
tic events and the associated elastic response, that re- 
laxes the stresses over the system. These minimal in- 
gredients are sufficient to generate a complex rheolog- 
ical behavior 15|. Approximating the rearrangements 



by spherical inclusions (Eshelby problem), one expects a 
four-fold quadrupolar symmetry for the inhomogeneous 
part of the stress propagator G{r,9) = ^^cos(46'). To 
simplify the model further we consider an incompress- 
ible medium and assume that the microscopic geometry 
of the plastic events is the same as the the one of the 
macroscopic shear, which permits a scalar description of 
the stress. The values of this local stress CTm are encoded 
on a square lattice of size N = L'^, given in units of a^, 
the typical size of a plastic event (e.g. several grains in 
case of granular material). 

The deterministic part of the stress dynamics reads, 
using dimensionless quantities 



(5) 
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where the time t is measured in units of an elementary 
relaxation time t, the stress a in units of the local yield 
stress Uy and the strain rate 7 in units of the critical value 
7c = '^yl'Tl^^ M being the shear modulus. 7c indicates 
the change from a Newtonian to a non Newtonian flow 
behavior. is the discretized propagator for a ffiiite 

geometry with periodic boundary conditions [isj . We as- 
sume for the dynamics of the plastic part of the strain 
eP' a viscoelastic like relaxation of the material, weighted 
with a local state variable indicating whether a site 
is plastically active or not e^^ = nia{t)aia{t)/2. The 
stochastic dynamics of the activity is ruled by the fol- 
lowing dynamics for the plastic state variable 
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In simulations we chose Tpiast = Toiast = 1 in units of 
r. Note that convection effects on the stress and the ac- 
tivity have been neglected to keep the model as simple 
as possible. The lack of convection leads basically to an 
additional symmetry for the spatial arrangements of the 
plastic events. This is expected to affect the shape of 
the heterogeneities, but will not alter the general rela- 
tion between dynamical heterogeneities and the diffusion 
dynamics. 

Dynamical heterogeneities. - To obtain the scaling of 
the cooperative volume Vc entering equation [31 we study 
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the dynamical correlations within the model described 
above. Under the influence of the external strain rate 
7, the system reaches a non equilibrium steady state, 
where two time correlations of all quantities in the sys- 
tem become time-translation invariant. The time needed 
to achieve this steady state increases as the inverse of 
the strain rate. A first indication for growing dynami- 
cal heterogeneities with decreasing strain rate is the be- 
havior of the two time correlation function of the local 
stress fluctuations pia{t) = <Jia — crd, where Ud is the 
steady state stress averaged over time and space. The 
correlation function reads C{t) — {pia{0)pia{t)/ p{0)'^)c, 
where the bar indicates a spatial average and the brack- 
ets a configuration average. In the limit of low strain 
rate, a data collapse is obtained when C is plotted as a 
function of strain (see Fig. Oi). The relaxation time is 
therefore inversely proportional to strain rate, although 
a deviation from this behavior is observed as the strain 
rate increases. A hint to the existence of a complex, co- 
operative dynamics lies in the non exponential behavior 
of C, which can be interpreted as resulting from the co- 
existence of regions with fast and slow dynamics. Such 
dynamical heterogeneities are revealed by the study of a 
four point correlation function, 

J/3 

- {Pia{0)Pia{t))c{Pjf:){0)PjfAt))c] ■ (6) 

The integral of G4 over space yields the variance of the 

two time correlation function, Xiit) ~ N{C2{t) — C'^(t)) 

2 

withC2(t) = {{pia{Q)pia{t)/ p{QY) )c- X4 geucrically dis- 
plays a maximum at a time t* corresponding to the time 
of the largest heterogeneity of the system, and in general 
is the typical decay time for C. Moreover, X4(<*) can be 
interpreted as the number of sites involved in events at 
the cooperativity maximum, and is therefore identified 
with the number Vc introduced in the first part. 

As expected from the study of C, t* scales with the 
inverse of the applied strain rate. A more interesting 
feature is the growth of Xi{i*) as 7 0, which clearly 
indicates the increasingly heterogeneous behavior of the 
system. A more detailed study of G4 shows that this 
function is strongly anisotropic, a result of the special 
form of the stress propagator (see Fig. Integrating 
G4 over space reveals a fractal dimension of the cooper- 
ative domains close to 3/2, which is consistent with the 
finite size scaling displayed in Fig. [2ji. 

To obtain this scaling we plot the height of the peak 
as a function of 7 for different system sizes in the inset of 
Fig. [2ji. Clearly we observe a strong size dependence. In 
an infinite system we expect that the fractal dynamical 
cooperation length, traditionally denoted as ^4(7), would 
grow as a negative power of 7. However, in a finite sys- 
tem, this length will saturate when reaching the limit of 
the sample. This effect leads to the plateau region in the 
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FIG. 2: (color online). Dynamical heterogeneities in the 
stress field. When not otherwise specified, the system 
size is A'^ — 2^®. Top left: Normalized two-time auto- 
correlation function of the stress as a function of strain 
(7 = 10"^10~^■^10~^10"^■^10"''). The inset shows the 
same quantity as a function of time. Top right: Time evo- 
lution of the dynamical susceptibility (same values for 7). 
Bottom left: Spatial shape of the normalized four point stress 



correlation function Ga for a strain rate 7 



10" 



at time 



t* where the dynamical susceptibiUty is maximal. Bottom 
right: Finite size scaling plot for the maxima of the dynam- 
ical susceptibilities. Shown are the maxima as a function of 
the inverse strain rate, both normalized by the system size (+ 
iV = 2^^ X AT = 2", * = 2^*^). The dotted line has slope 
one and the inset shows the raw data. 



Xi curve with value xo • For a fixed system size we reach 
this plateau below a critical strain rate 7^. It is then 
natural to write the following finite-size scaling 



XA{i.L)^L^'^f({Ui)lL) 



3/2 



(7) 



with a scaling function fix), that is linear for small val- 
ues of X and saturating for large x. This hypothesis is 
perfectly confirmed in Fig.[2Jl, with ^4 ^7"^. We expect 
the specific value for the fractal dimension for the geom- 
etry of the heterogeneities to depend on some specific 
details of the model. Here, cooperative events are free to 
spread in the directions of the two main axes due to the 
additional symmetry introduced through the absence of 
convection in the model system. This results into fractal 
structures with dimension larger than unity, rather than 
one-dimensional slip lines. We believe this behavior to be 
representative of three dimensional systems where an ad- 
ditional direction is introduced, in agreement with earlier 
results on 3D molecular dynamics simulations [l8| . 
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FIG. 3: (color online) Left: scaling of the difltusion coefficient 
D with strain rate and system size {* N = 2^ , + N — 2^^\ ■ 
N = 2^^). The inset shows T) versus strain 7 for 7 = 10^'^. 
Right: Master curves for diffusion coefficient ( x ) and dynam- 
ical susceptibility (+), normalized by their plateau value as 
a function of strain rate in units of the corresponding critical 
strain rate, indicating the onset of finite size effects. The inset 
shows the unsealed master curves. 



Diffusion ~ Within the present model, we introduee 
diffusion of tracers by considering the non-affine motion 
originating from the long range displacement fields in- 
duced by plastic events. We then associate with each 
(mesoscopic) event a corresponding continuous displace- 
ment field (as introduced in the former scaling argument) 
and use it to define the mean square displacement of the 
tracer particles. Nonintcr acting tracer particles arc as- 
signed to every lattice site, and their fictitious trajectory 
is built progressively by adding up contributions from all 
plastic events. In this way, the mean square displace- 
ment (Ar^) can be obtained and a diffusion coefficient 
T> = 72? with V = (Ar^)/2A7 can be extracted with a 
good accuracy. The mean square displacement as a func- 
tion of strain is ballistic for small strains, and becomes 
diffusive after a transient regime (A7 « 1), with a well 
defined slope D, see inset in Fig.[3K. This quantity varies 
strongly with strain rate (for strain rates 7 < 7c) and sys- 
tem size (for all 7), see Fig. For a given system size 
we find power law behavior for the scaling with strain 
rate D ^ •j^^-^ down to a critical rate •jo below which 
a plateau is reached. This plateau value Dq scales with 
system size as Dq ~ _L~^'^. To summarize, we find the 
same scaling behavior as for the dynamical susceptibility, 
see Fig.[3]D, with a crossover length scale for diffusion £d 
proportional to 7^^. 

Therefore we validate the scaling prediction in Eq. |4j 
-D ^ X4i foi' large values of Lj. Only the onset on the 
finite size effects occur for different parameter values, 
7£) w 37^ (see inset of Fig. (Sb). This indicates that the 
intrinsic length scales for the heterogeneities, ^4 and for 
the diffusion, io, are not identical but only proportional 
with ^4 < £d- 

Discussion. - Our analysis shows how the mechanism 
of deformation in amorphous materials at low tempera- 
ture, described by elastic interactions between isolated 
shear transformation events, results in nontrivial cooper- 



ative effects that make a strong contribution to passive 
tracer diffusion at low strain rates. This connection is 
embodied in a scaling relation between the diffusion co- 
efficient and the cooperative volume for plastic hetero- 
geneities. We stress that the diffusion considered here is 
associated with the long range part of the elastic defor- 
mation induced by plastic events, thus does not involve 
the more direct one, associated with the motion of parti- 
cles at the core of the transformation zones. The former 
is expected to be dominant at small strain rate. The 
peculiarity of this diffusion is a strong system size de- 
pendence, that reveals the non local aspect of the system 
dynamics. 

In the literature, several observation of strain rate and 
system size dependent diffusion constants have been re- 
ported, both in experiments [l2l . [l9| and simulations 
[id 1201 1 . Our results suggest that in addition to strain 
rate effects, the measured diffusion may be strongly af- 
fected by system size effects. It would therefore be desir- 
able to perform experimental finite-size scaling analyses 
of tracer diffusion in flowing amorphous materials, us- 
ing e.g. confined microfiuidic systems, which provide key 
information about the rheological heterogeneities. 
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